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Abstract 

Increasing requirements in the sensitivity of interferometric measurements is a common feature of 
several research fields, from gravitational wave detection to quantum optics. This motivates refined 
studies of high reflectivity mirrors and of noise sources that are tightly related to their structure. 
In this work we present an experimental characterization of photothermal and thermo-refractive 
effects in high reflectivity mirrors, i.e., of the variations in the position of their effective reflection 
plane due to weak residual power absorption. The measurements are performed by modulating 
the impinging power in the range 10 Hz 100 kHz. The experimental results are compared with 
an expressly derived theoretical model in order to fully understand the phenomena and exploit 
them to extract useful effective thermo-mechanical parameters of the coating. The measurements 
are extended at cryogenic temperature, where most high sensitivity experiments are performed (or 
planned in future versions) and where characterizations of dielectric film coatings are still poor. 



1 



INTRODUCTION 



Extremely sensitive interferometric measurements are implied in several frontier experi- 
ments. For instance, large baseline gravitational wave detectors presently surpass the limit 
of 10^^^ Hz~^/^ in relative displacement spectral noise jl|. An irnprovement of more 



than one order of magnitude is planned for their advanced versions and a sensitivity 
of 10~^^ m/A/Hz is envisaged in the kHz range for possible future massive, cryogenic de- 
tectors [2^. On the small scale, very low- noise background is also obtained in small scale 
quantum opto-mechanics experiments jsl- 

With such extreme requirements, critical limits are imposed by several noise sources 
associated with the interferometer optical components {4], and a crucial role is played in 
particular by the high-reflection coating. For instance, Brownian thermal noise of the di- 
electric layers is the most stringent foreseen limiting factor in the ~kHz range for the next 
generation of gravitational wave detectors. Another important phenomenon is the pho- 
tothermal dynamic effect: an (even small) fraction of the light impinging on a mirror is 
absorbed and heats it, causing a local temperature change and, through thermal expansion, 
a deformation of the mirror surface. A fundamental limit to the interferometer sensitivity is 
therefore produced by the shot noise of the absorbed radiation, and such limit can be worse 
in the case of excess of laser-power fluctuations. This effect, together with other dissipa- 
tive phenomena such as photoelastic and Brownian noise, has been studied theoretically by 
Braginsky et al. |4,] and Cerdonio et al. jsl in the case of a Gaussian beam impinging on a 
half- infinite homogeneous mirror. 

In spite of being an annoying noise source, photothermal dynamic displacement can 
be exploited, using controlled laser power variations, for measuring the thermo-mechanical 
parameters of the mirror. This possibility is particularly important to investigate the char- 
acteristics of the coating. Thin film parameters are generally different from those of the 
bulk. Moreover, they depend on the film morphology, and finally on the techniques used for 
deposition. A photothermal analysis gives therefore useful information for understanding 
the general behavior of the coating and selecting the optimal fabrication procedure. At this 
purpose, it is important to analyze in particular the low-losses coating and substrates of the 
kind used in interferometers. 

Dynamic photothermal effects were observed in a high-Finesse Fabry-Perot cavity, prob- 



ably for the first time, by An et al.\^. They describe an hysteretic behavior in the cavity 
transmission when the laser is rapidly scanned, and use it to estimate the mirror absorption. 
A more detailed spectral analysis exploiting an amplitude modulated laser source is reported 
by De Rosa et al. [7] , who verify the theoretical prediction of Refs. [4] and [5] concerning the 
photothermal expansion of the mirror substrate. 

The study is extended in Ref . 8j , where on one hand the dependence of the photothermal 
response on the beam waist, foreseen in Ref. [sl, is verified; on the other side, a deviation 
from the above theoretical expression is observed at low frequencies (< 10^^ Hz) due to the 
limited mirror size, and at high frequencies (> 10^ Hz) due to the coating. The effect of 
the coating is also shown by Black et a/.jsl, who apply a gold film on the dielectric coated 
mirror to enhance the absorption at the expenses of the cavity Finesse. Their dielectric film 
is composed of alternate layers of Si02 and Ti02, and their results are in agreement with a 
simple heuristic model approximating the coating as an homogeneous layer. 

In view of their possible use in the future gravitational wave detectors, it is important 
to check the photothermal effect also on substrates of different materials (particularly inter- 
esting are sapphire and silicon) and at cryogenic temperature. Such conditions, explored by 
the pioneering work of the TAMA group 10|, would greatly improve the sensitivity (at the 
expenses of evident technical complications). 

In the present work we describe a detailed study of photothermal dynamic effect in a high- 
Finesse cavity, where the mirror coating is very close to those conceived for gravitational 
wave detectors. Our attention is indeed particularly focused on the behavior of the dielectric 
films, but we also check the effect on a silicon substrate. Our investigation is performed 
both at room and at cryogenic temperature, and a theoretical model is developed in order 
to extract the most reliable information from the experimental data. 

MODEL 

The photothermal expansion is ruled by the thermal diffusion equation and the equation 
of elasticity. The mirror surface displacement is then read with a weighted average that, in 
the case of a Gaussian beam, can be written as 
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where Uz{r,t) is the surface normal displacement, r are the transverse coordinates of the 
mirror surface, vq is the beam radius measured in correspondence of a power drop by 1/e 
from its central value. 

The solution of the problem in the case of an half-infinite homogeneous mirror, heated 
and interrogated by a Gaussian beam partially absorbed on the mirror surface, is calcu- 
lated by Braginsky, Gorodetsky and Vyatchanin (BGV). A more conclusive form of the 
photothermal noise spectral density is given by Cerdonio et al. [5] in the framework of an 
extended analysis of thermoelastic noise sources. Their expression can be transformed to 
give the response to a modulated absorbed power 5Pabsi that in the Fourier space reads 

«H = -^!^G('^)^p„„M (2) 
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{5X{uj) is the Fourier transformed of u{t)) where a, a, A are respectively the thermal expan- 
sion coefficient, the Poisson ratio, and the thermal conductivity of the mirror; ujq = X/pCr^, 
where p is the density and C the specific heat capacity; the shape of the response is 

G(SJ) = - du jjo ■ (3) 

The requirements for the validity of the above expression (e.g., a short mean free-path of 



the phonons, and a large sound velocity) are discussed in details in Refs. 

At high frequencies {u ^ uq, the so-called adiabatic limit), G{Q) tends to This 
regime occurs when the thermal diffusion length in one modulation period is much smaller 
than the beam radius. At even higher frequencies, when this length is comparable with 
the coating thickness, the approximation of homogeneous mirror fails and a more complete 
model is necessary. 

For an immediate understanding of the physics involved, it is useful to consider a one- 
dimensional model, starting from the thermal conduction equation that reads 

dT d'^T 

where T is the temperature, the mirror occupies the space with z > 0, and Q is the absorbed 
power. In the frequency domain, the above equation can be written as 

k'-^)f = Q/X (5) 



where we have introduced a complex wave number 

k = 



iujpC 
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(6) 



The inverse of the absolute value of k gives the penetration depth of the thermal wave. 
In the case of an absorbed intensity labs on the surface at 2; = 0, the boundary condition 

is 

-A^(0) = J,,, (7) 

and, in a multi-layer medium, at each interface Zi the continuity conditions for temperature 
and heat flux are 



(8) 



T{z)l^=T{z)l, 

The general solution of the homogeneous equation associated with the ([5]) is 

T = Ae^' + Be-''' . (9) 

In the case of a mono-layer with thickness d and associated wave number /c^, lying on a 
substrate with wave number fc^, the temperature is given by Eq. ([9]) on the layer and by 
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(10) 



on the substrate, where the coefficients A, B, and Bg are derived from the boundary con- 
ditions and (IE]). The temperature on the layer can also be written in a compact form 
as 



rfnf ^ _ Lbs sinh[fce((j - z)] +ncosh[kc{d - z)] , , 
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where 
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and the subscripts s and c distinguish respectively the thermal parameters of the substrate 
and the layer. The pre-factor determining the temperature in the substrate is 

_ ^abs 1 

Xgkg cosh[A;c'^] + 'R- smh.[kcd] 
We remark that Eq. ( fT3l) differs from the solution in the absence of the layer (i.e., with 
= 0) by a filter function 
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cosh [/ccfi] -|- 7^sinh[fcc'i] 
5 



(14) 



which describes the modification of the heat fiux caused by its transmission through the 
mono-layer. 

The filter function can be generally defined as 

df 



(15) 



where in general the heating power can be absorbed both on the surface and within the 
medium: 



■Af(0) 

dz 



Qdz 



As suggested by BGV, the mirror surface deformation can be assumed to be 



5u, = - I a{z)T{z)dz. 
'o 



(16) 



(17) 



In this way, the one-dimensional solution differs from the adiabatic approximation only for 
the lacking term (1 + (x), if one uses for labs the peak intensity 6Pabs/ (tTq). 

It is useful to notice that integrating Eq. ([5]) along an homogeneous layer and using 
Eqs. ( !T5|) and ( !T6|) we obtain 



k^X I T = I. 
'o 
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This expression can be used to write the surface displacement in a compact form, even with- 
out performing the direct integration of the temperature along the layer. Indeed, calculating 
the mirror thermal expansion according to Eq. ( IT7|) and using Eq. (fT8|) and the exponentially 
decreasing temperature in the substrate, with the pre-factor of Eq. ( IT3l) . we obtain 
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(19) 



The above equation is correct in the one-dimensional approximation, that can be relaxed 
assuming a Gaussian beam with radius tq ^ d, an absorbed intensity labs = 5 Pabs / {T^fD , 
and correcting the expansion coefficient of the substrate with the factor (1 + cr). Observing 
that k"^ = (i/rD^u/uo) and defining Uc,s similarly to coq, we can write 



6X{u) 



TcXr iuj ttAo ioo 



5Pabs{(^) ■ 



(20) 



On the right side of the above expression, the first term is due to the expansion of the layer, 
while in the second term we see the expansion of the substrate caused by the filtered thermal 



flux. Besides the filter, the second term coincides with the adiabatic limit of Eq. (|2]), as in 
the case of an homogeneous mirror. A smooth transition to the homogeneous, 3-dimensional 
model can therefore be simply obtained by replacing the factor {ug/iuj) with the function 
G{uj/(jJs) defined in equation (E]). 

A more rigorous expression for 5X{ijj) is calculated in the Appendix A using a 3- 
dimensional formalism for both the thermal conduction and the elastic equations. The 
calculation considers an half-infinite mirror with a superimposed layer, thus neglecting the 
finite mirror dimensions. The result is summarized in Eqs. fl7m74p . 

One further thermo-optical phenomenon that must be considered is the thermo-refractive 
effect: heating changes the optical length of the coating layers, thus modifying the phase 
of the reflected field and therefore the position of the effective reflection plane with respect 
to the physical position of the mirror surface. We remark that an increased optical length 
gives a longer effective cavity (while the thermal expansion brings to a shorter physical 
cavity), therefore the thermo-refractive effect has in general an opposite sign with respect to 
photothermal expansion. The attention to the thermo-refractive effect has been drawn by 



Braginsky et al. [13|, who identify it as source of thermo-elastic noise. A correct expression 
for the thermo-refractive coefficient in the case of our interest (i.e., a coating starting from 
the high index layer, without cover layer) is reported in Ref. [1^. Since the effect is mainly 
due to the first few layers, we can adopt the simplified expression valid for infinite number 
of layers and using the temperature at the coating surface: 

= X,ac(3efff{0) (21) 

where the thermo-refractive coefficient is 

^'^^^ = i^l^ ^^^^ 



with Bh,l = ^ + aH,LnH,L uM- In the framework of a one-dimensional model, the 



temperature T(0) can be taken from the expression ffTTj) using z = 0. As for the thermal 
expansion, the contribution to the signal 6Xtr{uj) is calculated using the spatial average 
given by Eq. ([T]). A more accurate calculation in a three-dimensional framework is reported 
in Appendix A, and the result is summarized in Eq. ( 1771) . 

So far we have only considered surface absorption. In fact, the light is absorbed in the 
coating, where the electromagnetic field forms a standing wave, with nodes/antinodes in 



correspondence of the interfaces between layers and a rapidly decreasing amplitude. More- 
over, the absorption is in general different in the two materials composing the coating. The 
exact solution of the field within the coating can be calculated using a characteristic matrix 



formalism 12|]. For a high- reflectivity coating (with many layers), the average intensity in 
the first layer is reduced with respect to the intracavity value by a factor very close to the 
refractive index nn of the material (in our case, Ta205). In each following layer, the intensity 
is further reduced with respect to the previous one by a factor ni/nn, where rii is the index 
of the second material (for us, Si02). If we consider that each layer is long (0.25 \vac/nH,L), 
where Xyac is the vacuum optical wavelength for which the coating is designed, we see that 
the absorption in the first period is Pcav ■ 0.25 X^ad^/^H + where luilh) is the 

extinction length in the H (L) material, and P^av is the intracavity power. In each following 
period, the absorption is reduced by a factor ini^jnu)^ ^ therefore the total absorption in one 



1 1 \ 1 



mirror is Pahs — AbsPcav where 

Abs = 0.25 A.,, + ■ (23) 

For our refractive index values of uh = 2.05 and ul = 1.46, about 74% of the absorption 
occurs in the first two periods. In brief, we observe that a) in the two mirrors the absorption 
is the same, since the kind of substrate and the number of layers are not critical parameters 
to determine the intensity decay within the coating; b) the two material composing the 
coating have the same weight in evaluating the total absorption, i.e., the absorption in each 
material is just determined by its extinction length, independently from its thickness; c) 
for our parameters, the approximation of surface absorption is well satisfied since, even at 
the highest experimental frequencies, the thermal penetration depth corresponds to several 
periods. 

EXPERIMENT 

Apparatus and procedures 



The optica 
fused sihca 



cavity is made with a OFC (Oxygen Free Copper) spacer, a Corning 7980 
concave mirror (50 mm radius) as input mirror and a flat, silicon rear mirror. 
Mirrors diameter and thickness are respectively 12.7 mm and 6 mm. The cavity optical 
length is 10.7 mm, the beam waist on the flat mirror is 0.083 mm and the corresponding 
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beam size on the input surface is 0.094 mm. The input mirror has a coating composed 
of alternate quarter- wave layers of Ta205 and Si02 (their physical thickness is respectively 
130 nm and 180 nm). The overall number of layers is 29, giving a transmission of about 
100 ppm and a total coating thickness of about 4.5 fim. The rear mirror has a similar 
coating, with 34 layers (5.3 /im) and a lower transmission (~ 10 ppm). 

The cavity spacer is a thick box, with just two 10 mm diameter holes for the optical 
access, and few smaller holes for evacuating. It is fixed in thermal contact with the cold 
finger inside a continuous flow ^He cryostat which is evacuated down to 5 • 10^'* Pa. 

Temperature sensors are placed on the cold finger and on the cavity spacer. To assure 
an efficient cooling and temperature stability of the cavity mirrors, they are mounted inside 
the OFC box. A spring keep them in contact with the box through their rear, fiat surface. 
The standard thermal shield of the cryostat has four, 25 mm diameter access holes. We 
have further reduced the apertures to few mm, with aluminum foil. As a consequence, the 
mirrors 'see' the room temperature through couples of holes (on the spacer and the shield), 
under a small solid angle. 

The experimental setup for the measurement of the photothermal response is sketched 
in Fig. ([T]). The light source is a cw Nd:YAG laser operating at 1064 nm. After a 40 dB 
optical isolator, the laser radiation is split into two beams. On the first one (probe beam), 
a resonant electro-optic modulator (EOMl) provides phase modulation at 13.3 MHz with 
a depth of about 1 rad used for the Pound- Drever- Hall (PDH) detection scheme. The 
probe beam is then frequency shifted by means of two acousto-optic modulators (AOM). We 
use respectively the +1 and -1 diffracted orders of the AOMs, so that the total frequency 
displacement corresponds to the difference of the AOM frequencies and can be tuned by 
several MHz around zero. The intensity of the second beam (pump beam) is controlled by 
an electro-optic modulator (E0M2) followed by a polarizing beam splitter. Both beams are 
sent to the second part of the apparatus by means of single-mode, polarization maintaining 
optical fibers. Part of the pump is then detected by a photodiode (PD2) for monitoring its 
amplitude modulation, while the probe is sent through a second optical isolator (O.I. 1). 
The two beams are overlapped with orthogonal polarizations in a polarizing beam-splitter 
and sent to the optical cavity. The refiected probe beam, on its back path, is deviated by 
the input polarizer of the 0.1.1 and collected by a photodiode (PDl) for the PDH detection. 
The pump, slightly detuned with respect to the cavity resonance, is used as excitation beam 
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in the thermo-optical response measurements. Such frequency sphtting between pump and 
probe allows to eliminate any spurious interference and reduce the cross-talk between them 
in the photo-detection. A quarter-wave plate before the cryostat allows to compensate 
cavity birefringence observed at low temperature and probably due to stress on the mirrors 
generated by differential contractions. The laser is weakly locked on the cavity using the 
PDH signal while the pump beam is amplitude modulated at different frequencies. The 
extraction of the photothermal response is obtained from the same PDH signal, that is 
acquired by a lock-in amplifier which discriminates the signal component synchronous with 
the modulation. The depth and phase of the intracavity power modulation is calculated 
from the signal of PD2, which is also acquired and de- modulated by the lock-in amplifier, 
by taking into account the mode matching and the cavity coupling factor. The latter is 
inferred from the depth of the dip in the reflected beam when scanning on the resonance, 
and it is consistent with the independent measurements of cavity Finesse and input mirror 
transmission. 

The PDH signal is calibrated in a separate run, stopping the pump beam and modulating 
the frequency of the laser through its internal piezoelectric crystal, with a depth smaller than 
the Fabry- Perot cavity linewidth. The usual lock-in acquisition and phase-sensitive detection 
is applied to the PDH signal. The frequency modulation rate had been previously measured 
by observing the sidebands at 13.3 MHz with a second, high-Finesse optical cavity. The 
calibration procedure, performed varying the modulation frequency, allows also to correct 
the data for the servo loop. The locking bandwidth is of the order of 1 kHz. 

Typical values of the optical parameters are: a mode-matching of 80%; a Finesse of 25000 
(linewidth ~ 475 kHz); an intracavity optical power Pcav equal to 7000 times the input pump 
power Pin ~ 2 mW. The input probe beam power is approximately 0.2 mW. 



Experimental results and data analysis 



The photothermal response to intensity modulation of the intracavity field has been in- 
vestigated both at room and cryogenic temperatures in the 10 Hz 100 kHz frequency 
range. The measured signal is reported in terms of changes of the optical path divided 
by the amplitude of the modulation in the intracavity power, giving a response in m/W. 
Previous experiments [7|, [Sj have shown that, at low frequencies (below ~100 Hz), the pho- 
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tothermal expansion is determined by the substrate contribution. In agreement with the 
theory of homogeneous mirrors, the curve can be roughly divided into two regions, sepa- 
rated by the characteristic angular frequency Ug which depends on the thermal properties 
of the mirrors material. Well above Ug-, the surface displacement is expected to decrease 
as the inverse of the frequency while, well below Ug-, it tends to a constant value which 
determines the absorption coefficient. On the other hand, at higher frequencies (above 
~100 Hz), the coating contribution to the total thermal expansion, in particular through 
the filter function, becomes relevant and the observed behavior differs substantially from the 
substrate 1// frequency dependence. At even higher frequencies, also the thermo-refractive 
effect becomes important (in the following, we will use the term 'photothermal response' for 
the overall effect, including the thermo-refractive component). An accurate data analysis 
through Eq. f lTHj) and the related expressions in Appendix A allows us to to measure the 
effective thermal properties of the Ta205/Si02 dielectric layers. Moreover, assuming that 
the thermal parameters of silica films are not significantly different from those of bulk silica, 
we will estimate in the next section the properties of Ta205 thin film. On this material, to 
our knowledge, just few, room temperature measurements have been performed. 

The amplitude and phase of the photothermal response at room temperature are shown 
in Fig. (|2]). The experimental data (amplitude and phase) are globally fitted through the 
complex response function ( |78l) considering both the silica and silicon substrate contributions 
and the coating. The fitting function is shown in Fig. (j2]) (red traces), and Fig. (I2ti) also 
shows the amplitudes of some specific contributions. In the fitting procedure, the substrates 
thermal properties (thermal conductivity, specific heat capacity and expansion coefficient) 
are fixed to values taken from the literature and reported in Table I. The coating parameters 
derived from the fit are reported in Table II (first column). 

The mechanical unbalance parameters 71 and 72 (see Eq. (1701) ) are both close to unit and 
depend on the Young modulus {E) and the Poisson factor (a) of substrate and coating. We 
have taken the values of E and a from the literature (they are reported in Tables I and II), 
including the experimental uncertainty for the coating, and a spread due to the dependence 
on the direction with respect to the crystal axes for Si, since we have no information on their 
orientation in the mirror. The calculated parameters are 71 = 1.36±0.05 and 72 = 1.53±0.12 
for the silica mirror; 71 = 0.97^0.99±0.02 and 72 = 0.58-^0.80±0.06 for the silicon mirror. 
The uncertainty in the mechanical parameters gives a small spread in the fitted parameters 
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(below the statistical error given by the fit), that we have included in their quoted uncertainty 
in Table II. 

The absorption coefficients of the two mirrors are kept as different parameters even if 
the two coatings are very similar (their slightly different thickness is not expected to have a 
great inffuence on the absorption), since they have been provided in different batches and we 
have no precise information on the post-coating annealing procedure that can be different 
for the two mirrors. 

The coating parameters (pcCc), etc, and /3e// are determined by the fit. Its results 
are reported in Table II together with their uncertainties derived from the 68 % confidence 
intervals, slightly increased to account for the spread in the used values of Young modulus 
and Poisson factor. 

At frequencies below 100 Hz the silica substrate contribution (whose characteristic fre- 
quency is ~ 27r ■ 27 Hz) is dominant while at higher frequencies the photothermal 
response is dominated by the sihcon substrate (wf * ~ 2n ■ 4200 Hz) and the coating. The 
importance of the coating is not restricted to its direct contribution to the photothermal ex- 
pansion (dash-dotted curves), but it is evident above all in the filter function that strongly 
affects the total response. This is put into evidence in Fig. ([2]d), where the phase data 
are also compared with a curve obtained form Eq. ([2]), thus neglecting the filter function 
contribution (dashed curve). In this case, the disagreement between theory and experiment 
is relevant even at low frequencies and would lead to unreliable estimations of both the 
substrate and coating parameters. The phase of the response is sensitive also to the thermo- 
refractive effect, which seems negligible if one just consider its amplitude. This is underlined 
again in Fig. ([2]d), where the dotted curve neglects 6Xtr in Eq. fl78l) . By considering the 
complete model (red curve), instead, the agreement is excellent in the whole frequency range. 

We observe that the coating on silicon mirror displays a lower absorption. The absorption 
coefficients that we have obtained are a bit higher than previously reported |^ , yet remaining 
within few ppm. We remark however than the quoted uncertainty is reliable for comparing 
the two mirror, but the estimate of the overall absorption is conditioned by the evaluation 
of the intracavity power that is hardly very accurate. 

The measurements have been repeated at low temperature, with the cavity spacer at 
4.5 K. The results are shown in Fig. ([3]), together with the fitting functions, and the de- 
rived parameters are reported in Table II (second column). Around this temperature, the 
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photothermal response is dominated by the fused silica substrate at all frequencies. 

We remark that, at low frequencies, the phase curves reported in Fig. ([2]d) and Fig. ([3]d) 
clearly show that the photothermal response at room and cryogenic temperature are 180° out 
of phase. This is due to the sign change in the fused silica expansion coefficient, which be- 
comes negative below ~ 170 K. The effect of the coating filter function on the photothermal 
response phase is even more important than at room temperature, as evidenced in Fig. ([3b). 

For what concerns the silicon substrate, on one hand the silicon thermal diffusion length 
becomes larger than the mirror size for frequencies below few kHz, thus invalidating our 
theoretical model (the temperature of the silicon substrate could be considered as homoge- 
neous, instantaneously determined by the heat absorption and the conduction through the 
copper spacer); on the other hand, its thermal expansion coefficient becomes so low that it 
can be completely neglected. Some figures justifying this assumption are reported in Table I. 

Concerning the thermal parameters of the silica substrate, we remark that expansion 
coefficient and capacitance have a strong dependence on temperature. Since the real tem- 
perature of the mirrors is also determined by their thermal contact with the copper spacer, 
we have kept it as free parameters in the fit, using the known temperature dependence of 
Csi02 and asi02 (reported in the caption of Table I). The derived temperature of the mir- 
rors is indeed 10.9 ± 0.3 K, giving for bulk sihca asi02 = -0.267 ± 0.012 ■ 10"^ K'^ and 
Csi02 = 5.2 ± 0.5 J/kg K. 

The mechanical parameters (71, 72, a) are kept at their room temperature values (their 
weak temperature dependence is negligible in our work). We have verified that their uncer- 
tainty has negligible effects on all the fitted parameters. 

The contribution of the coating on silicon is taken into account, imposing 7?. = in 
the relevant equations as justified by the silicon thermal parameters reported in Table I. 
The ratio between the absorption coefficients of the two mirrors is fixed at the value found 
at room temperature (i.e., Ahssi = QAA Ahssi02)- The ~ 50% uncertainty in this ratio 
adds a small error to some fitted parameters (well below the corresponding statistical errors 
form the fit), that is included in the uncertainties quoted in Table II. The absorption in 
the silica mirror results to be lower than at room temperature. However, we remark again 
that the quoted uncertainty in the absorption coefficients only reflects the statistical error, 
while absolute intracavity power calibrations (giving additional systematic uncertainty) are 
different between runs at room and cryogenic temperature due to the necessary re-alignment. 
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The thermo-refractive effect at cryogenic temperature is too small to allow a reliable 
independent estimate of both ac and /3e// from our data. Yet, it cannot be completely 
neglected. Therefore, we have left etc as free parameter and used the expression /3ejj of 
Eq. fl22]) as follows: for the Bi^ parameter (concerning silica) we have used the silica bulk 
values of a and dn/dT reported in Table I; the expansion coefficient of tantala is derived 
from Eq. (126|) using the free parameter ac and the bulk expansion coefficient of silica; we 
have varied the thermo-optic coefficient dn/dT of tantala between and twice the silica 
value. The spread in this coefficient has negligible effect in all the fitted parameters except 
for where the spread is comparable to the statistical uncertainty of the fit. Since our 
assumption on the range of the tantala thermo-optic coefficient is rather arbitrary, we have 
quoted separately its effect on Uc in Table II to allow a reliable evaluation of our results 
once an independent measurement of dn/dT for tantala is provided. 



Discussion 



As discussed before, the fitting procedure allows the estimate of the coating effective 
thermal conductivity Ac, the product pcCc-, the averaged thermal expansion coefficient ac and 
the thermo-refractive coefficient /3e//- The values of these effective parameters are the main 
result of our work, and they can be directly used in most evaluations of noise sources and, 
in general, physical phenomena related to the high reflectivity coating in interferometers. In 
this section, we further discuss our results with the aim of extracting the parameters of the 
materials composing the coating, and comparing with existing literature. 

The models used in our data analysis consider indeed an homogeneous coating layer. In a 
complete description of the film behavior, the real parameters of the layers would enter with 
several different combinations that should be individually averaged over the coating. How- 
ever, to extract some estimates of the layers parameters, we adopt the simplified description 

I \ 

suggested, e.g., in Ref. [ill, that takes volume averages such as 

+ (24) 



Ac dp \Xh X^ , 

pCc = Y (^HPH Ch + dLPL Cl) (25) 
dp 

andu + aidi .„„x 
ttc = -. ■ (26) 

dp 
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dp = dH + di (27) 

where the subscripts 'if' and 'L' refer to the two consecutive layers composing one period 
of the coating. Eqs. f l2Hl271) and fl22l) allow us to extract the parameters of the materials 
composing the film from our experimental results. 



We start our discussion from the room temperature data. Ref. 201] presents a review of 
conductivity measurements in Si02 and Ta205 films, for different thickness and deposition 
techniques. The range of results is very wide, extending between 0.017 and 0.73 W/mK 
for silica, and between 0.026 and 15 W/mK for Ta205. If we assume, in spite of that, for 
the Si02 layers the bulk parameter values reported in Table I, from Eq. (!24|) we obtain for 
Ta205 a thermal conductivity of ATa2O5=0.4 W/mK. We stress that such estimate should 
be taken cum grano salis, since conductivity in thin films can strongly depend on the film 
structure; the most meaningful parameter remains Ac. 

Concerning the heat capacity, it can be estimated using the bulk parameters for silica, 
the measured density of 7200 kg/m^ for the Ta205 film 18|, and the bulk tantala capacity 
of 306 J/kg K |l9|. We calculate from Eq. ([25]) a value of (p^Cc) ^ 1.9 ■ 10^ J/m^ K. The 
result of our fit is (pcCc) = 1.4 ± 0.3 ■ 10^ J/m^ K, that is compatible (within two standard 
errors), while a bit lower, with the previous estimate. 



The thermal expansion coefficient of Ta205 can be extracted through Eq. (12611 . using 
again the bulk data for silica. We obtain ara2O^=4:.S±0.7 ■ 10~® K~^ to be compared with 



23| measured in previous 



Tom the analysis of data 



22 



231 than to the one 



the values of -4.4 10"^ K"^ [21], 2.42 10"^ K^^ [22] and 5 10^^ K"^ 
experiments. Ref. 22| also report a value of 3.6 10^^ K^^ obtained 
in the literature. Our value is closer to those reported in Refs. 
obtained by Inci 21| which, apart the unexpected negative sign, is more than a factor of 10 
larger in absolute value. As discussed in Ref. |23|, the variability between these independent 
measurements could be due to the fact that the properties of tantalum pentoxide layers 
might depend on the procedure of the layer deposition. 



An accurate measurements of the re 
from 30 to 310 K is reported in Ref. 



:r active index of bulk silica at temperature varying 
24j |. From its data, we infer (at X^ac = 1064nm) 



dn/dT = 8.57 ■ 10"^ K'^ at 300 K, and dn/dT ~ (8 + 3.3T) ■ 10"^ K^^ at low temperature. 
The room temperature value is compatible with dn/dT = 1.2 ■ 10~^ K~^ measured at 
632.8 nm for a sputtered SiO^ film in Ref. 25|]. Using Eq. ( 122|) and the bulk data for 
silica, we obtain for the tantalum pentoxide dn/dT + na = 27 ± 5 ■ 10~^ K^^, and from our 
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value of aTa205 we calculate {dn/dT)Ta205 = 17 ± 7 ■ 10 ^ K ^. Previous measurements of 
{dn/dT)Ta205, performed at 632.8 nm b y a group in Taiwan, gives 2.3-10"^ [25| and, 
later, 3.64-10-'' K'^ and 7.89-10-'' K" 



value of 14-10 ^ K ^ is mentioned in 



26| (depending on the thermal stress). A further 



15| as a private communication. 



At low temperature, the literature on the coating properties is very poor. By using, 
as before, the expressions (!24|) - (l271) together with the bulk properties of Table I for silica, 
we deduce the following data for the film of tantalum pentoxide: the expansion coefficient 
is 0^0205 =0.58±0. 04 ± 0.05 - 10^^ K^^ (the first quoted uncertainty is due to the spread 
in the assumed thermo-optic parameter of tantala between and 2 x dn/ dTsi02i with 
dn/dTsi02 = 0.44 - 10^^ at 10.9 K; the second quoted uncertainty is the statistical error 
from the fit); the thermal capacitance is {pC)Ta205 = 2.2 ± 0.7 - 10^ J/K m^ and, using a 
density of 7200 kg/m^, we calculate 6*^0205 = 3.1 ± 1.0 J/kg K; the thermal conductivity is 
XTa20s = 5 ± 3 - 10-2 W/mK. 



CONCLUSIONS 

We have measured the variation of the effective reflecting surface in a high reflectivity 
mirror, caused by a modulation in the residual absorbed power. Such complex response is 
due to photothermal and thermo-refractive effects in the mirror substrate and coating, and 
our measurements allow a characterization of these effects both at room and at cryogenic 
temperature. An accurate knowledge of such phenomena is particularly relevant in experi- 
ments aiming to the highest sensitivity in displacement, position and force measurements, 
for instance in the flelds of gravitational wave detection and in quantum optics. Most of 
the planned future experiments of this kind will address cryogenics to reduce Brownian 
noise, and part of the present work is indeed the focused on cryogenic temperature, where 
measurements of the kind here presented were still missing. 

Our experimental results are well described using a theoretical model that we have ex- 
pressly developed, and the comparison allows us to extract meaningful effective thermo- 
mechanical parameters of the high reflectivity coating. This is the most important aspect of 
our work, since a) the overall coating behavior is not obviously deduced from the properties 
of the composing materials (in particular concerning thermal expansion and conductivity); 
b) the parameters that we have measured can be used to calculate weak but important 

16 



noise contributions such as thermodynamic fluctuations or, of course, photothermal and 
thermo-refractive noise. 

We under hne that, while specific experiments on the different oxides in the bulk or thin- 
film form can give more accurate results concerning the material properties, our setup has the 
added value of characterizing the mirrors (and coating) as they are used in high sensitivity 
experiments, thus describing their real overall behavior. Moreover, from the effective coating 
parameters we could still estimate some properties of the dielectric films composing the 
coating. In particular, literature on previous measurements were rather poor for tantalum 
pentoxide at room temperature, and, in our knowledge, completely lacking at cryogenic 
temperature. 
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APPENDIX A 

In this Appendix we calculate the temperature profile and the elastic deformation of an 
half- infinite mirror extended in the half-space z > 0, composed of a layer of thickness d an 
an homogeneous substrate, with surface absorption of a Gaussian beam. We will consider an 
isotropic substrate, neglecting the directional variation of the Si properties. The calculation 
is performed using Fourier transforms for the temporal and transverse spatial variables, 
while keeping the dependence on the longitudinal variable z. 

Thermal equations 

We start by the study of the temperature distribution. We define 




(28) 



and 




(29) 
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pC dT . _ 
A ot 



where r = (x, y, z),r± = {x, y), and kj_ = {kx, ky). The equation for the thermal propagation 

(30) 

can be written in the time-transformed Fourier space as 

k'^f-V'^f = Q (31) 



where k is defined in Eq. (I6l), and in the space- and time-transformed Fourier space as 



[e + ki)Q- 



(32) 



The absorbed intensity, in the time-transformed Fourier space, is written as 



/(u;,rj_l 



(33) 



and, in the complete Fourier description where X is the spatial Fourier transformed of /, it 
becomes 



X(a;, k^) = exp 



k^: r,^ 



(34) 



The three-dimensional model is now formally equivalent to the one-dimensional model 
described in Section „ The equations fl32l) is equal to the homogeneous equation associated 
to Eq. ([5]), and the boundary conditions ([7]) and ([8]) are still valid, with the following list of 
substitutions: 

f ^ e (35) 



labs — ^ X 

k^^iP^co + kl = k'' + kl 
A 



(36) 
(37) 



The discussion on the temperature distribution in the one-dimensional model remains valid 
too. We can therefore use Eqs. fl MT^ . with the replacements fl35ti37p and where kg and kc 
contain the respective thermal parameters. For the sake of clarity, we repeat here explicitly 
the solutions, in the layer (from Eq. flTTj) ): 



e 



exp (^—^-^J ^^^^ 


^kl + kl{d- z) 


+ TZ cosh 


+ kl{d- z) 




cosh 


^kl + kld 


-|- TZ sinh 


^kl + kld 





5PaU^) (38) 
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and in the substrate (from Eqs. fll0|13l) ): 

exp 

e = — ^ 



exp 


-^kl + kl {d- 






cosh 


^/kl + kld 


+ 71 sinh 


^kl + kld 



(39) 



where the reflection coefficient is now 



7^ 



K^/W+ki 



(40) 



Elastic equations 



We now consider the problem of the mirror thermal expansion. The elastic deformation 
u is ruled by the static elastic equation 



^V(V-u) - i^^V A(VAu) = aVT 
1 + cr 2(1 + a) 

The relevant components of the stress tensor are defined as 

E 



(41) 



1 - 2a 



V ■ u + 



1 + a 
E 



l + a dz 



2(1 + a) \dx dz 
E ( duz ^ du, 



(42) 

(43) 
(44) 



2(1 + a) \dy ' dz 
where E is the Young modulus and cr is the Poisson ratio. 

Following BGV, we search a particular solution of Eq. f l4T|) in the form of a gradient of 
a potential (p. Inserting V0 in Eq. (HTl) . we see that the potential (p obeys the Poisson's 
equation 



l + g 
l-a 



a T 



We then search a complete solution of the kind 

u = V + V0 



(45) 



(46) 



that is inserted in Eq. (HTi) to show that v satisfies the associated homogeneous equation, 
that we recast in the form 



V(V ■ v) + (1 - 2a)VV = 
19 



(47) 



As in the study of the thermal diffusion, we define the time- and transverse spatial 
Fourier transformed variables kj_, z) of 0; Ui^u, k±, z) and Vi{uj, k_L, z) of the Cartesian 
components of u and v, respectively, with i = x,y,z. The equation (H7|l . in its Cartesian 
components, can be written in the Fourier space, giving 

d / dV \ ( 92 \ 

— [tKV^ + ^k,V, + ^ j + (1 - 2 a) [--kl + ^ J K = (48) 

from the z— component, 

\ik^V^ + ikyVy + ^ j + (1 - 2 a) f -fcl + ^ j \4 = (49) 

from the x— component, and the equivalent of Eq. ( H9l) with exchange of x and y for the 
y— component. 

The axial symmetry of the problem implies that Vx/k^ = Vy/ky. We can therefore define 

Vr = = t^Vy (50) 



that, replaced in Eq. fl48l) . gives 
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^^^ + 2(l-^)(-^l + ^jK = (51) 

and replaced in either Eq. ( H9|) or its equivalent with x y (the two equations transforms 
into the same expression) gives 

-kJ^ + 2{l-a)(-kl + ^\Vr = Q . (52) 



dz ' ^ ' \ ^ ' d z"^ y 
We remark that the condition V^/kx = Vy/ky simplifies very much the calculation. How- 
ever, one can also avoid to impose it at the beginning, and verify later that it results as a 
consequence of the boundary conditions. 

Eqs. f l5T]) and f l52l) form a system of homogeneous, second order differential equations with 
constant coefficients, in the variable z, for the vector of functions {Vi,V^r}- The solution is 
a linear combinations of the four vectors {e^^^ , —e''-^^}; {e~^^^,e~^^^}; {k±ze^^^,{4a — 
3 — k±z) e^^^}; {k± ze~'^^^, (4cr — 3 + fc^ z), e~^^^}. It is useful to express the solution in 
different forms for the layer (m Q < z < d) and in the substrate (for z > d), where we 
skip the diverging solutions with positive real part of the exponent. Therefore, we write the 
solutions as 

Vz = Vc cosh{k±z) + Vs smh{k±z) + Vcck± z cosh{k±z) + Vgs k±_ z smh.{k^z) (53) 
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Vr = —Vc sinh(fcj_z) — Vg cosh.{k±z) + Vcc [(4ac — 3) cosh(fc^z) — k±z smh.{k±z)^ 

+Vss [(4(Tc - 3) smh{k±z) - k±z cosh{k±z)] (54) 

in the layer and 

V, = VaC-^^^''-^^ + Vbk^ {z - d) e-^^(^-'^) (55) 

Vr = Va e- + Vb [{4as -3) + k±{z- d)] e" (56) 

in the substrate. Here dc and cr^ are the Poisson factors respectively in the layer and in the 
substrate, and Vc, Vg, Vcc, Ks, K, H are functions of u and k±. 

For what concerns the particular solution V0, instead of solving Eq. ( H5|l . we can compare 
it with Eq. fl3T]) and directly write (in the complete Fourier space) 

1 + CT 1 



$ = a ■ 



e 



(57) 



1 - a F 

The definitions of the stress tensor given in Eqs. fl42ti44p can also be written in the 
complete Fourier space. Naming dij the space- and time- Fourier transformed of aij, using 
Eqs. (146|) and (!57|) . we obtain 

E 



I -2a 



a ^ , 1 — a dVz 1 — 2a A , n ^ 
k± Vr + ^ + a 



a 



1 + a dz 



E 



a 



2(1 + a) 



1 dVr 



+ 2 



1 + a 



a iu pC 
A dQ 



■ a ■ 



(58) 



(59) 



kj_ dz 1 — 0" iu pC dz 
Once written a complete solution of the elastic equation, we now consider the boundary 

conditions. At z = Q the null stress condition is 



azz\ n = 







=0 ^^y|z=0 

The continuity conditions for the stress at the interface between layer and substrate are 



(60) 
(61) 



I z=d 



a-. 



I z=d~ 



=d+ 



The continuity of the displacement at the interface is given by the relation Ui\^_^_ 
that can be transformed, using Eqs. fH6!) and (1571) . into 



Vz + a, 



1 + ac 1 dQ 
1 ~ ac kl dz 



z=d~ 



V, + a, 



1 + a, 1 dQ 
1 — (Jo /c? dz 



(62) 
(63) 

M2=<i+' 

(64) 



=d+ 
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Vr — ki a. 



l+(Tr 1 



1 - a, k^. 







z=d 



Vr - k± as iJi-^ -L e 

1 - ex., kj 



(65) 



z=d+ 



Replacing the solutions ( I38ll39|) and fl53tl56p into the Eqs. (I60tl65p we obtain a system of 6 
linear equations for the 6 variables Vc, Vg, Vcc, Vgs, Va, H, that can be solved using annoying 
but straightforward algebra. As we will see, for our purpose the meaningful function is just 



Vr.. 



Measured thermal expansion 



The average displacement measured by the laser beam is given in Eq. ([T]). By using the 
inverse Fourier transform for Uz{r±,t) we can write 



u{t) = / due"^'6X{u 



e 
Tirl 



d^ki 



2 = 



(66) 



and performing the integration over and over the angular coordinate of , we can express 
the time-Fourier transformed of such displacement as 

dfc 



5X{lo) 



27r 



/c_L exp 



z=0 



(67) 



Using Eqs. ( I46l) . (1531) and (|571) the function Uz in the integrand of Eq. ( 1671) is transformed 
into 



(68) 



2=0 



The expression of | (and, in particular, of V^) is completely and explicitly determined 
by the algebra described above. It can be shown that Uz \ can be written in the following 
form 



Uz 



2ar 



2=0 ^2 



+ /c_L sinh(/c^(i) G _, — cosh(A;_L(i) 



2 = 



dz 
(96 

+ 72(fc±) ( /^±0|^^o - kx_ cosh(fc^ci) e|^^^ + sinh(fc^(i) — 



z=d 



+ -^73(fc. 



dz 



+ k^Q\ 



=d+ 



(69) 



The particular combinations of the temperature functions that appear in Eq. fl69l) assure 
the correct frequency dependence for w — )■ 0, compensating the divergence given by l/k^g. 
Therefore, they must be kept as they are. On the other hand, the functions 71 (/c±), 72 (/^±) 
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and 73(fc±) are rather complex for a practical use. It is useful to give their approximations 
at the lowest order in k\d, that read 



7i 

72 



7i(0) = - 



72(0) 
73(0) 



21 -a, 

1- (JcEs 

1 + as 



1 + (1 - 2(T, 



(70) 



These three parameters are close to unit, to which they reduce if as = Cc = and Ec = 
Es- When inserting f^z|^_Q in Eq. fl67|) . the approximated expressions of the 7's can be 
used, provided that d « ro/2, thanks to the fast decaying Gaussian function of (fc^ro/2). 



Replacing 9 from Eqs. fl38ti39p in Eq. ( |69l) . inserting the resulting Uz\.^^q in Eq. (!67jl . 
and after some simplifications, we can write the expected photothermal response to the 
modulated absorbed power as the sum of two terms, equivalent in the simplified model to 
the expansions of coating and substrate: 



- 72 



7i ^cosh(^d/ro) + 7^|- sinh(^ci/ro) -cosh(^crfAo) - TZsinh^^cd/ro] 
^ 7^cosh(^ d/ro) + ^ sinh(^ d/ro) — Tlcosh{^c d/ro) — smh.{^cd/rQ] 



and 



(71) 



(72) 



where the adimensional integration variable is ^ = k±rQ, and we have defined to simplify 
the notation 

'kl^ + P, = .A— +e • (73) 



The filter function F, derived from Eq. ( fT4l) by replacing kc — )■ a/^^ + can be written 

as 

F = . (74) 

cosh [(.cd/rol +71 sinh [C,c d / Tq] 

In the limit d = (no coating layer) the filter F obviously assumes a unit value, and 6Xs 
coincides with the results of Refs Q] and [sl (reported here in Eq. ((21)). In other words, while 
not evident at a first glance, it can be shown that 



1 



1 - 



(75) 
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where G{Q) is given in Eq. ([3]). Concerning the expansion of the coating, it can be shown 
that 6Xc in the high frequency hmit reduces to the one-dimensional expression (first term 
on the right side of Eq. (I2CT]) ). if one set 71 = 72 = 0. 

One further effect to be considered is the thermo- refractive displacement. From Eq. (pTj) . 
using Eq. fl67j) for the average displacement, we can write it as 







6Xtr = KacPeff / ^ exp ( ) (76) 



and replacing 6 from Eq. fl38|) 

Xvac/3eff r sinh(^cC?/ro) +7^ cosh(^crf/ro) , . 

27r Ac ro Jq 4c cosh (4c a/ro) + Tc smh(4c d/ro) 

Finally, the total response to the modulated absorbed power is 

6X = 6X, + 6Xs + 6Xtr . (78) 
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SiOs (T = 300 K) 


SiOs (low T) 


Si (T = 300 K) 


Si (T = 10 K) 


Cs (J/kg 


K) 


790 


Csi02 (T) 


710 


0.4 






Ref. [16] 


Ref. [27] 


Ref. [33] 


Ref. [35] 


A. (W/m 


K) 


1.30 


0.12 


150 


1600 






Ref. [16] 


Ref. [27, 28] 


Ref. [32, 3^ 


Ref. [32] 


as (K-M 




5.2 10"^ 




2.6 10-6 


4.8 10-1° 






Ref. [16] 


Ref. [29-31] 


Ref. [34] 


Ref. [36] 


Ps (kg/m' 




2201 
Ref. [16] 




2330 




E (GPa) 




72 

Ref. [16] 




130 [100] 

L J 

169 [110] 
188 [111] 
Ref. [37] 




a 




0.16 
Ref. [16] 




0.28 (100) 
0.21 (110) 
0.18 (111) 
Ref. [37] 








8.57 10-6 ( 
Ref. [24] 


;8 + 3.3 T) 10" 
Ref. [24] 


-8 





TABLE I. Substrate parameters used in the fitting procedure. The expressions for the low- 
temperature thermal capacitance and expansion coefficient of silica are respectively Csi02{T) = 
(0.64 + 0.54r - 0.02ir2)r3 • lO-^ and asio^ (T) = (2.6 - 5.2r + 2.9r2 - OmT^ + 0.044T'i) • 10"^ . 
The Poisson factors for Si are averaged over the respective planes. 
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T = 300 K low T 



PcCc (J/m^ K) 1.4±0.3 -10^ 1.6±0.3 -lO"^ 
Ac (W/m K) 0.7±0.3 7±3 -IQ-^ 
Uc (K-i) 2.3±0.3 -10-6 9±1.5±2 -lO"^ 

Peff (K~i) 4.4±0.6 -10-6 
Mssi02 (ppm) 9.0±0.5 3.8±0.2 
Ahssi (ppm) 4±2 0.44 ■KAhssi02 

E (GPa) 91±7 

Ref. [38j 
cr 0.20 

Ref. [38] 

T (K) 10.9±0.3 



TABLE II. Coating parameters in the fitting procedures. For Uc at low temperature, the first 
quoted uncertainty reflects the spread in the assumed thermo-optic coefficient of tantala. 
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FIG. 1. Scheme of the experimental apparatus. O.I.: optical isolator; AOM: acousto-optic modu- 
lator; EOM: electro-optic modulator; H: half-wave plate; Q: quarter-wave plate; PD: photodiode; 
PBS: polarizing beam-splitter; BS: beam-splitter. 



29 




FIG. 2. a) Amplitude of the photothermal response as a function of the modulation frequency at 
room temperature. The solid red line is the fit through the complex response function (j78p for 
both silica and silicon mirrors, keeping fixed the substrate parameters (see Table I). Other curves 
refer to the fused silica (green) and silicon (blue) mirrors. Dashed lines: substrate contributions; 
dash-dotted lines: coating expansion contributions; dotted lines: thermo-refractive contributions, 
b) Phase of the photothermal response and its corresponding fit. Blue curves are traced neglecting 
some effects, in order to show their relevance: the dashed line corresponds to a response of just 
homogeneous substrates (calculated from Eq. the dotted line considers just 5Xs and 6Xc in 
Eq. ()78p . thus neglecting the thermo-refractive effect. 
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FIG. 3. The same as in Fig. ([2]), at cryogenic temperature. 
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